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We study one-dimensional anisotropic XY-Heisenberg spin-^ chain with weak random fields /if Sf 
by means of Jordan-Wigner transformation to spinless Luttinger liquid with disorder and bosoniza- 
tion techniqne. First we investigate phase diagram of the system in terms of dimensionless disorder 
7 = /J^ ^ 1 and anisotropy parameter A = Jz/Jxy and find the range of these parameters 

where disorder is irrelevant in the infrared limit and spin-spin correlations are described by power 
laws. Then we use the diagram technique in terms of plasmon excitations to study low-temperature 
behavior of heat conductivity k and spin conductivity a in this power-law phase. The obtained 
Lorentz number L = nfaT differs from the value derived earlier by means of memory fnnction 
method. We argue also that in the studied region inelastic scattering is strong enough to suppress 
quantum interference in the low-temperature limit. 


I. INTRODUCTION 

One-dimensional disordered spin chain is an excellent 
example of strongly correlated quantum system that is 
well suited to study basic properties of such systems. In 
particular, studies of disordered spin chains become one 
of the major playgrounds in the field of Many Body Lo¬ 
calization (MBL) — . From the experimental viewpoint, 
quasi-one-dimensional antiferromagnets^— attract con¬ 
siderable attention due to very high thermal conduc¬ 
tance, that is believed to be related with integrability 
of the clean Heisenberg spin-i chai n^^d^ . It is known 
since seminal paper— that in ID competition between 
interaction and disorder may lead to delocalization and 
formation of a ground state that is nearly-free from the 
effects of disorder, see also ^^d^ . Numerical studiesii con¬ 
firm that qualitative conclusion. In order to provide de- 
localization, interaction should be sufficiently strong and 
attractive, so this problem bears some resemblance with 
a model of superconductor-insulator transition in higher¬ 
dimensional systemsiS. Looking from that perspective, it 
seems useful to develop a quantitative theory of the delo¬ 
calized phase of one-dimensional quantum system with a 
bare disorder that is ’’screened” by interactions. In par¬ 
ticular, it is important to study heat transport in such a 
system, that is expected to be dominated by the remains 
of the disorder potential. 

Here we will study the properties of anisotropic XXZ 
spin chain in a random transverse magnetic field, which 
is described by the Hamiltonian (we assume J > 0): 

H = ( 1 ) 
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By means of the Jordan-Wigner (JW) transformation 
the Hamiltonian ([T|) can be reduced to the Hamiltonian 
of interacting spinless fermions (here = cJjCn, — i): 

H = J ^ ^ h.C ^PnPn+1 ( 2 ) 

The anisotropy parameter A can be both positive and 


negative, which corresponds to the effective attraction or 
repulsion between JW fermions, respectively. 

We will consider random fields to be relatively small so 
that 7 = (/i^)/J^<Cl, and with zero average (h) = 0. 
Thus our system is, on average, symmetric with respect 
to z —z reflection, which translates into the particle- 
hole symmetry in terms of JW fermions. It ensures that 
in the quasiparticle spectrum £(k) only odd powers of k 
survive. 

The goal of this paper is to study low-temperature 
transport properties, spin and heat conductivity, in the 
range of parameters (A, 7 ) where T — 0 spin-spin corre¬ 
lations decay as a power law with a distance. The rest 
of the paper is organized as follows: in Sec.H we study 
the phase diagram by means of Renormalization Group 
approach formulated ii>i^. Sec.HI is devoted to formu¬ 
lation and application of the Keldysh approach to the 
transport properties of disordered Luttinger liquid model 
that is an appropriate low-energy approximation for the 
lattice fermion model (ED; in Sec.HI A, spin and heat 
conductivities (cr and k) are studied within the region 
i < A < 1 where disorder is irrelevant in the RG sense; 
next, in Sec.HI B we discuss specific behavior of a and k 
near the critical point A = i; the role of quantum inter¬ 
ference corrections and decoherence is discussed in Sec. 
HI C, and the role of spectrum nonlinearity is considered 
in Sec.HI D. Finally, we present our Conclusions in Sec. 
IV. 


II. LUTTINGER LIQUID DESCRIPTION AND 
PHASE DIAGRAM 

In the clean limit /i„ = 0 and in the region — I < A < I 
excitation spectrum of the interacting one-dimensional 
fermion system ([2) is gapless; then low-energy and long¬ 
distance properties of the system are known to be de¬ 
scribed by the Luttinger liquid (LL) modeli^. It al¬ 
lows to rewrite the Hamiltonian in terms of fermion 
density excitations — plasmons. LL model is formu¬ 
lated in terms of canonically conjugated plasmon fields 
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n( 2 /)] = i5{x — y); in the linear approximation for 
the quasiparticle spectrum, the Hamiltonian of the LL 
model reads 

=-^ J . (3) 

Here u is plasmon velocity and K is dimensionless Lut- 
tinger parameter; these parameters are determined, via 
the Bethe Ansatz solution for XXZ model, by the values 
of J and A, seei^ (p.l67): 


A = 


cos ■ 


2K' 


Ja sin( 7 r/ 2 Ar) 

T 1 - 1/2A: 


( 4 ) 


where a is the lattice constant. 

In our model (HI) disorder couples to the fermion den¬ 
sity pn', in the LL continuum limit it reads as p(x) = 
— cos{2kFX — 2(^). First and second terms in 

the above expression correspond to the slow {q ^ 0 ) and 
fast oscillating {q ^ 2A:f) parts. Thus there are two types 
of scattering of one-dimensional fermions by disorder: 
forward and backward. Forward scattering is irrelevant 
within the linear approximation for the spectrum, since 
the corresponding term in the LL Hamiltonian can be 
eliminated completely by the redefinition of phase ^(a;). 
Backward fermion scattering with momentum transfer 
q ~ 2kF is the only effect one should take into account 
then. Thus we need only q ~ 2kF part of original random 
potential, this part is described by the random Gaussian 
complex field ^{x) with {^{x)^*{y)) = DS{x — y) and 
D = a. Disorder contribution to the Hamiltonian 
reads as follows: 

Hd^s = J + r (5) 


Renormalization Group approach to disordered Lut- 
tinger liquid was formulated ini^. It is convenient to 
introduce dimensionless disorder parameter 


_ 2Da _ 8(1 - 1/2K)^ 

^ 7^2 7rsin2(7r/2A:)^ 

In terms of this parameter and logarithmic scaling pa¬ 
rameter ^ = In ^ (with a being running ultraviolet cut¬ 
off), RG equations reads as follows: 


du 

di 

dK 

di 

dg_ 

di 


-^9 



= {3-2K)g. 


( 7 ) 


These equations can be solved analytically exploiting 
its first integral I{K,g) = |(A-|-41nAr — y). This so¬ 
lution yields phase diagram shown in Fig. [TJ ’’Delo¬ 
calized” region lies, in the limit of very weak disorder, 
in the range ^ < A < 1. Upon increase of 7 , delocal¬ 
ized region shrinks and eventually disappears already at 
7 Ri 0.1. Everywhere in the delocalized phase effective 
disorder y(^) decreases with Actually the derivation 


of the RG equations © as it was performed ini^ is valid 
quantitatively in the vicinity of the point IF = 3/2 only, 
where disorder-induced corrections to the parameter K 
are logarithmic; for large K these equations can be used 
for qualitative analysis only. Note that the drop of the 
critical disorder 7 near the point A = 1 is trivially re¬ 
lated to the decrease of the effective Luttinger velocity 
M, see Eqs. 0 . 

Phase diagram obtained by the analysis of RG equa¬ 
tions can be compared with the numeric phase diagram 
from RefAZ (its boundary is shown by the dashed line in 
Fig. H]). According to these numerical data, the delocal¬ 
ized region covers much smaller part of the phase dia¬ 
gram than the RG calculations predict. We expect that 
the major source of this discrepancy is due to inapplica¬ 
bility of the RG equations ([7]) at large K values. Another 
reason could be related with the effects of spectrum non¬ 
linearity that becomes important close to A = 1. On the 
other hand, near the point K = 3/2 numerical dataiii 
suggest delocalization at the values of 7 which are above 
our critical line; we believe that this discrepancy comes 
from limited accuracy of the numerical data, due to finite 
size effects which becomes most prominent at very weak 
disorder. 



FIG. 1. Approximate phase diagram found from the RG 
calculations for the Luttinger liquid model with linear spec¬ 
trum. Dashed line corresponds to phase boundary obtained 
in Ref.— . 


Equal-time spin-spin correlation function 
(S'+(0)S'“(a;)) decays as a power law at A > 4, as 
one can read of Refp^^ where two-loop RG calculation 
was performed. At smaller A < 4 renormalized disorder 
parameter y(^) grows with and one expects exponen¬ 
tial decay of {S^(Q)S~{x)) at x > Lc, where correlation 
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length Lc cx seei^. 

Note that for the case of XY model with random 
transverse fields (i.e. A = 0) exponential decay of 
(S'+(0)S'“(a:)) follows directly from single-particle local¬ 
ization in ID, as proven rigorously in Ref;^. However, 
at general A < | the relation between growth of effective 
disorder upon RG and Anderson localization is far from 
being obvious, since the RG calculation — does not con¬ 
tain any multiple-impurity interference effects, se o^°'^^ . 

Below we will focus on delocalized phase i < A < 1, 
that corresponds to the range of | < K < oo, where 
renormalized disorder constant g is small, and one can 
obtain transport properties using perturbation theory for 
bosonic LL model with renormalized parameters. 


III. TRANSPORT PROPERTIES 


Here we proceed from the Hamiltonian description de¬ 
fined by Eas. (l3l5p to the Keldysh action for the LL model 
with disorder. Total Keldysh action Stot consists of triv¬ 
ial free boson part So and disorder-related part coming 
directly from Eg. ([5]) : 

Sd^s = T^I (8) 

27ra Jc 

We integrate exp(iS'tot) over random Gaussian field ^(x) 
and perform Keldysh rotation introducing classical (pci = 
^{(/>+ -I- (j)-) and quantum pq = (j)+ — </)_ fields compo¬ 
nents, arriving finally at the effective disorder action 


iD f 

Sdts = -5-^ / dtidt 2 dx cos 2 { 4 >ici - 02d) sin^i, sin(/)2q 

(9) 

In order to obtain self-energy for retarded Green func¬ 
tion in the lowest order over Sdis we consider first order 
correction to it, which reads 

iSG^]^\y - y') =J dtidt2 j dxi2 x 

X {4>cl(l>'qCOs{2{(t)lcl - 4>2ci))s.\n(l)iqS\n4>2q)^ (10) 

Here we used notation y = (y,t) and pi = p{y,ti). Per¬ 
forming Wick’s contraction, one finds two diagrams (see 
Eig. [^, from which we extract retarded bosonic self¬ 
energy Efl;(w); the corresponding analytical expression 
reads 

4D 

E«(a;) = -^/ dt(l-e“‘)x 
0, Jo 

xe2*(°^(‘)-GK(o))sin2Gfl(t) (11) 


where bare retarded and Keldysh components of the 
Green function are as follows 




ttuK 


(uj + iO)^ — 
ttK , 

g’-rIuj) = coth^(G®(a;) - G^°^(w)) 


G^°\t,x) = - \x\) 


( 12 ) 

(13) 

(14) 



FIG. 2. Lowest order diagrams for the retarded self en¬ 
ergy T.r{ljj). Dashed lines correspond to disorder average 
(C(®)G(y)) ^nd wavy lines correspond to averaging of cosine 
or sine of boson fields 


Inverse Eourier transformation of the Keldysh compo¬ 
nent G^^ (w) to the real space-time, G^^ {t,x), is infrared- 
divergent; it is sufficient to use the difference G^^\t, x) — 
G^^(0,0) which is finite: 



2o2 




G^°\t,x)-G^°\0,0) 

. , irix + ut) . , irix — ut) 

Sinn-Sinn- 

up up 


(15) 


At low temperatures T J two different types of con¬ 
tributions to the disorder-induced self-energy can be sep¬ 
arated: virtual transitions with T uj < J and real 
(dissipative) transitions with oj < T. First contribution 
lead to logarithmic renormalization of the model param¬ 
eters yielding RG equations © described above; second 
contribution yields dissipative behavior of corresponding 
self-energy E/j = —iuj/ utiKt. Direct calculation yields 
the following expression for momentum relaxation rate: 


1 _ 2DKT'^{K) /2TraTy^~^ 

t{T) ~ u T{2K) V u ) 


(16) 


According to Ea. dTBl) . product Tr oc T3-2ic (diverges as 
T —>■ 0 in the delocalized phase. Full Green function then 
reads as follows: 


GR{q,u}) 


ttuK 

LtP- — -f fw/r 


(17) 


A. Spin and heat conductivities 


To obtain transport properties, one can apply Kubo 
formulas. Expressions for spin and energy currents can 
be derived from corresponding continuity equation dtPa + 
Vja = 0 (index a corresponds to either spin or energy), 
and using classical equations of motion. For the Hamilto¬ 
nian of the form H = J dxpE{Pix),'Vp{x),Il{x)), equa¬ 
tions of motion reads as follows: 


dtp 


dpE 

9n ’ 


dtli = 


dpE ^ dpE 
dp diyp)) 


(18) 


so energy density obeys the following continuity equation: 


dtPE 


O X I ^P^ O V7X I ^P^ O TT V7 

-dtp+i::;^dtVp+^:^dtIl = V 


dp 


dVp 


dW 


/ dpE dpE\ 
\dVp SH ) 
(19) 
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FIG. 3. Loop diagram for “thermal susceptibility” xe- Ther¬ 
mal current vertices act as the following combinations of 
derivatives: (9tV' -I- dfV), with derivatives dt, V acting on 
one <j) field in the vertex and d[, V' acting on another one. 


and similarly for spin density. Considering total Hamil¬ 
tonian consisting of two contributions and , we 
arrive at the following expressions for currents: 

1 u 

js = Je = - —dt4>V(t>. ( 20 ) 

TT TTA 

We emphasize that Eas. (l20t provide exact (within Lut- 
tinger liquid approximation) expressions for both spin 
and thermal currents. Surprisingly, in the LL approx¬ 
imation the energy current does not contain any terms 
related to the presence of backscattering. In Appendix 1X1 
we provide a detailed derivation of the energy current, 
starting from the lattice fermion model o, and show 
that backscattering does produce additional terms for the 
energy current, but these terms vanish in the continuous 
LL limit, when a —>■ 0 at some fixed value of the product 
Ja. 

Spin transport is governed by the single-plasmon Green 
function, while for energy transport we need to calcu¬ 
late correlation function of four (j) fields. Applying Kubo 
formula for spin conductivity, we reproduce Drude-like 
result of Refsi ^^i^^ . 

q = 0) = —-—-— (21) 

TT^ TT 1 — lUIT 

valid at a; <C T. 

Thermal conductivity k is expressed in terms of 
so-called “thermal susceptibility” XE{q,i^) as k,{u}) = 
~ X£;(0;0))- Introducing short notation 
X = {x,t), and q = {q,uj), expression for thermal suscep¬ 
tibility in real space reads as follows: 


Comparison between Eas. (l21l24l) provides us with the 
value of the Lorentz number 


L = 


K 


SK 


(25) 


2 

which matches its standard Fermi liquid value Lee = ^ 
for AT = 1. Note that our result (|2^ differs from one 
obtained in 22 by means of memory function formalism. 
We believe that this discrepancy is due to limitations of 
the memory function formalismwhich is based on ex¬ 
trapolation from the large-w region to the static limit. 
Indeed, frequency-dependent thermal conductance k{uj) 
depends on two different frequency scales, T and 1/r; 
according to Ea. IfTO . in the region K > "ij2 one always 
has T ^ l/r(T) in the low-temperature limit. In or¬ 
der to obtain static thermal conductivity, one should be 
able to compute k{uj) at wr <C I, whereas memory func¬ 
tion method is based upon the calculation of the high- 
frequency limit k{uj > 1/T) and further extrapolation to 
zero frequency. We believe that the presence of two para¬ 
metrically different frequency scales l/r and T makes 
such an extrapolation unreliable. 

The above calculation leading to Eas. d^ and (1^ 
should be performed, in general, with the parameters 
{g,K) renormalized (due to RG equations ([7])) down to 
the temperature scale ^t- If bare parameters {g, K) are in 
the bulk of delocalized phase (not too close to the transi¬ 
tion line) one can neglect renormalization of K and u due 
to disorder, leading to the results for spin and thermal 
conductivities which depend on the scale via scatter¬ 
ing time T only, see (HU). Then the result is given by Eqs. 
(HH), dH with bare parameters. 

Near the transition line one should take renormaliza¬ 
tion of all the parameters simultaneously. Below we will 
see how it affects physical properties of the system. 


B. Vicinity of the point A = 4 

Expanding first integral of system (0 by a: - or, 
equivalently, A — i, one obtains: 


X£;(xi - X 2 ) = i- 


.2K2 




and x(q) is the Fourier transform of this expression. 

In the dc limit w —^ 0 one finds k = —i/3 
Applying Wick’s theorem, one finds: 


a ;=0 


XB(q') = 


2ujq — 


X ^^(q + %)G^(q — ^) 


(23) 


Calculating it in the a; —?► 0 limit with “dressed” Green 
functions, we arrive at: 



^sinh^^ " ^ 


(24) 


I-Ic 



(26) 


The equality I = Ic yields the phase boundary of the 
delocalized state in the form (A — i) = ^ 7 . 

Solution of the equations ® can be expressed in terms 
of vicinity to transition line a = — 1^ <C 1 : 


ru(^) = uexp(|A:(^) - 1) 

<1 A:(^) = acothQ!(^-f ^ 0 ) (27) 

9 sinh^ 

where ^0 depends on initial values of parameters. Consid¬ 
ering temperature to be low enough (so ^ I'Coli l/o^); 
one obtains low-temperature behavior of renormalized 
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parameters gi^r) — a^exp(— 2 a^) ~ a^(T/J)^“ and 
K{^t) — I — a = const. 

Now we repeat the above calculations leading to 
nonzero ImE(w) and obtain Drude-type formulae with 
corrected power-law exponent a: 

CT ~ a-^a{T/J)-^-^°‘ K ~ Ja-^a{T/J)-^°‘ (28) 

The Lorentz number is still given by Ea. (l25l) once renor¬ 
malization K — 7 > K{^t) is taken into account. Modifica¬ 
tions of K and a are negligible if 5 <C (ill — 


C. Smallness of the interference corrections. 

Our result for the heat conductance, Eq. (HI, was ob¬ 
tained within Drude-type approximation. Since our sys¬ 
tem is one-dimensional, some care should be exercised to 
check if the effects of quantum interference and Anderson 
localization could affect that result. To begin with, it is 
useful to employ the result of Refi^ where the same is¬ 
sue was considered for disordered Luttinger liquid with a 
weak interaction, jAT — 1| <C 1. Namely, it was found in2£ 
that interference corrections are negligible at sufficiently 
high temperatures T > t~^{T){K — 1)“^. We are work¬ 
ing at K > 3/2 and the corresponding condition is just 
T ^ 1/t(T) which is always fulfilled at low temperatures 
according to Eq. dUD. 

To estimate interference corrections more accurately, 
we examine expression for “thermal susceptibility” to 
higher order in Sdis adding impurity lines connecting up¬ 
per and lower Green functions drown in Eig.3. Eirst or¬ 
der correction (with single impurity line) vanishes at zero 
external momentum due to gradient structure of energy 
current vertex. First non-trivial corrections are due to 
diagrams shown in Fig. 01 the corresponding analytical 
expressions yield: 


SI 1 

(5x£:(x-x) = - 


2 V TT^a^ 


dtidt2dt3dt4 X 


X y (i2;i2da;34^9i(()ci(x)V(/ci(x) X 

X (5i(/ci(x')V(/q(x') -f at(/q(x')V(/ci(x')) X 

X cos(2(/ici - 2(j)2ci) sinc/i, s\n(j) 2 q x 


X cos( 2 (/) 3 ci - 2(j)4ci) sin (/a, sin (j) 4 g ) (29) 


Grayed box correspond to sine and cosine average and 
consists of infinite number of boson propagators con¬ 
necting all the points. Generally speaking, such box 
depends on all the ingoing energies and momentums. 
However, direct calculation shows that it contains fac¬ 
tors l/sinh^^ 7 rT(t, - tj), which 

impose effective constraint for time differences: any such 
diagram is very small unless the condition — < 1/T 

is fulfilled. On the other hand, typical time scale for the 




FIG. 4. Non-trivial corrections to “thermal susceptibility”. 
Dashed lines correspond to same impurity, and grayed area 
correspond to average of cosine and sine of (f) fields. See main 
text for the analytical expressions. 



FIG. 5. Effective form of the diagram for the vertex correc¬ 
tion to thermal susceptibility, valid in the leading order of 
expansion over 1/Tt{T) ^ 1 


dressed “external” (w.r.t. to the ’’grey area”) propaga¬ 
tors is t{T) l/T; therefore, up to the leading order in 
1 /Tt(T) -C 1 one can try to shrink all four space-time 
impurity points in Fig. 0] into single one (see Fig. 0]). 
However, calculation of the remaining integrals result in 
a zero result, due to vector structure of the current ver¬ 
tex. Therefore, nonzero vertex corrections appear in the 
next order in 1/Tt{T) <C 1 only, and are small at low T 
in the whole “delocalized” phase K > 3/2. 


D. Spectrum nonlinearity effects 

At the Heisenberg isotropic point A = 1 in the clean 
system the spectrum of excitations is quadratic and sys¬ 
tem is no longer described by Luttinger liquid model. 
In the vicinity to this point plasmon velocity vanishes 
as M = Jav'(l-A)/ 2 ; since dimensionless disorder 
strength g depend on velocity u and interaction param¬ 
eter AT, see Eq.(|ni), this narrows the region where per¬ 
turbation theory in powers of small g is applicable to 
(h^) IJ^ < (1 - A)3/2. 

However spectrum nonlinearity effects at finite temper¬ 
atures might become relevant long before A = 1 critical 
point. Let us make some estimates. Due to particle- 
hole symmetry, only odd powers in quasiparticle spec¬ 
trum survive; first non-vanishing contribution to disper¬ 
sion relation will be 5e ~ ^{ka)^. At finite tempera¬ 
tures this yield new energy scale Se ^ T{Ta/u)'^; such 
energy scale should be compared with scattering ratio 
l/r ^ ug/a{T/JY^~^. Therefore we conclude that spec¬ 
trum nonlinearity will be important and should be taken 
into account when {T/ > (h?) j,P. 

For A" < 5/2 it leads to the threshold for the tempera¬ 
ture, above which nonlinearity effects are expected to be 
important, T* ^ ; on the contrary, 

at AT >5/2 nonlinearity is always important at low tem¬ 
peratures. In terms of the A parameter, the borderline 
at AT = 5/2 corresponds to A = cos 7 r /5 = (1 + •\/5)/4 ~ 
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IV. CONCLUSIONS 

We have analyzed spin and thermal conductance of 
XXZ spin chain with random-field disorder in the param¬ 
eter region where major source of disorder (backscatter- 
ing of Jordan-Wigner fermions) is suppressed by quan¬ 
tum fluctuations and irrelevant in the RG sense at low 
temperatures. Within the standard bosonization scheme 
the problem is reduced to the Luttinger liquid model with 
linear spectrum e{k) ~ uk and Luttinger interaction pa¬ 
rameter K in the range 3/2 < K < oo, which corresponds 
to 1/2 < A < 1 in terms of original anisotropy parameter 
A = JzjJ of the spin chain. We derive a phase bound¬ 
ary in terms of A and normalized disorder Ih?') / and 
compare it with the numerical result of Ref.f-. Then we 
use diagrammatic Drude-like calculation of thermal and 
spin conductivities and found Lorentz number, (1251) . dif¬ 
ferent from the previous result^^. We also argue that 
quantum interference (the effects beyond Drude approx¬ 
imation) is irrelevant at low temperatures due to strong 
enough inelastic scattering atl/2<A<l. 

These results were obtained neglecting forward scat¬ 
tering of Jordan-Wigner fermions by disorder, which is 
allowed as long as the approximation of LL model with 
linear spectrum is employed. However, this approxima¬ 
tion is not evidently correct everywhere in the delocalized 
phase. We estimated region where it might lead to quali¬ 
tatively different low-temperature behavior as A > cos ^. 
The effects of spectrum non-linearity will be considered 
in the separate publication. 

We have not studied the region A < i where localiza¬ 
tion due to disorder is expected; here it is very interest¬ 
ing to consider the close vicinity of the transition point, 
i — A <C 1 and (h?^ /<C 1 and search for the exis¬ 
tence of localization-delocalization threshold as function 
of excitation energy, like the one studied i n^^i^^ for the 
Bethe lattice model. 

This research was supported by the Russian Founda¬ 
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eral phase diagram, Fig.2) and by the Russian Science 
Foundation via grant # 14-42-00044 (all other results, 
presented in Eos. (1241251281) 1. 
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Appendix A: Derivation of the energy cnrrent 
starting from the lattice model. 

Eor a general Hamiltonian which is a sum of local on¬ 
site energy operators H = hn,n+i, with on-site ener¬ 
gies which satisfy continuity equation dthn,n+i +jE,n+i — 


jE,n = 0 with energy current jE,n = i 

For particular Hamiltonian (ED on-site energies has the 

form: 


^n,n+l 


= - J 


QX QX 


qy qy 






(Al) 

Substituting this expression into the expression for en¬ 
ergy current yields two contributions. First contribution 
is of kinetic nature, it does not contain disorder and can 
be written in a compact form as a determinant: 


■ (kin) 

3e,u 


= det 


qx 

^n-1 

qx 

‘^n 

qx 

•^n+l 

qy 

‘^n-l 

Sv 

qy 

•^n+l 


qz 

A5^+i 


(A2) 


Below we will focus only on the second term, which con¬ 
tains disorder. Corresponding expression in the original 
spin representation and in the Jordan-Wigner represen¬ 
tation reads as follows: 


.(dis) 

JE,n 


- s-_,st) = 


(A3) 


Next step is to take continuum limit by replacing lattice 
operators Cn with continuous field 4’{x = na) = Cnl\fa 
and replacing fields with continuous potential V(x = 
na) = hn- Corresponding expression for energy current 
density then reads as follows: 


= -i^V(x){'ip^x - a)4’{x) - 'il^\x)'ip{x - a)) 

(A4) 

In order to separate forward and backward scatter¬ 
ing, we introduce slowly varying in space left- and 
right-moving fermionic fields ipL,Rix) with = 

e''^^^ii}R[x)+e~'^^‘"^4’L[x). After splitting potential V{x) 
onto “forward-scattering” part "q^x) with Fourier har¬ 
monics g ^ 0 and “backward-scattering” part ^{x) with 
q ^ 2kF, one obtains contributions to energy current 
from forward- and backward-scattering processes: 

■U-s.) ^ _ a)i)R{x) - il)\{x - a)il)L{x)) + 

+ h.c. Ri Jaq{x) {i)'^p.{x)iiR{x) - %1 ;\{x)iI)l{,x)'^ (A5) 


■{h-s.) ^ _ a)ij}R{x) - V’l{(a:)V’L(a: - a)) -k 

-I- h.c. r: + V'|j(a:)VV’L) + h.c. (A6) 

One can see that backward-scattering contribution is of 
the next order in the small lattice constant a and indeed 
vanishes in the continuum limit, that is a —>■ 0 keeping 
u (X J a constant. 
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